Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN MfeK style file vl.4) 



Black— Hole Transients and the Eddington Limit 

A.R. King 1 

1 Astronomy Group, University of Leicester, Leicester LEI 7RH, U.K 



1 February 2008 



o 
o 
o 

(N 
X> 

IX, 



> 

00 
(N 
(N 
O 

o 
o 

a,; 

6 

in 



ABSTRACT 

I show that the Eddington limit implies a critical orbital period P cr it(BH) ~ 2 d 
beyond which black-hole LMXBs cannot appear as persistent systems. The unusual 
behaviour of GRO J1655-40 may result from its location close to this critical period. 



1 INTRODUCTION 



It is now well understood that the accretion discs in low- 
mass X-ray binaries (LMXBs) are strongly irradiated by the 
central X-rays, and that this has a decisive effect on their 
thermal stability (van Paradijs, 1996; King, Kolb & Burderi, 
1996). Irradiation stabilizes LMXB discs compared with the 
otherwise similar ones in cataclysmic variables (CVs) by re- 
moving their hydrogen ionization zones. In CVs this instabil- 
ity causes dwarf nova outbursts, and in LMXBs it produces 
transient outbursts rather than persistent accretion. The ir- 
radiation effect appears to be weaker if the accretor is a black 
hole rather than a neutron star, possibly because of the lack 
of a hard surface (King, Kolb & Szuszkiewicz, 1997). The 
result is that neutron-star LMXBs with short (~ hours) or- 
bital periods tend to be persistent, while similar black-hole 
binaries are largely transient. Both types of LMXBs must 
be transient at sufficiently long orbital periods, since a long 
period implies a large disc, so that a large X-ray luminos- 
ity would be needed to keep the disc edge ionized and thus 
suppress outbursts. We can write this stability requirement 



3 4/3 



(1) 



where M^rit i s the minimum central accretion rate required 
to keep the disc stable, is the outer disc radius, and P is 
the orbital period, and we have used Kepler's law. Thus for 
large P, M" r \ t will rise above any likely steady accretion rate, 
making long-period systems transient. This simple predic- 
tion (King, Frank, Kolb & Ritter, 1997) seems to be borne 
out by the available evidence. 



2 THE CRITICAL ACCRETION RATE 

The precise coefficient in (jl|) depends on uncertainties in the 
vertical disc structure (see the discussion in Dubus, Lasota, 
Hameury & Charles, 1999). Here I adopt the form derived 
by King, Kolb & Szuszkiewicz (1997). They argued that for 
a steady black-hole accretor, the central irradiating source 
is likely to be the inner disc rather than a solid spherical 
surface, as for a steady neutron-star accretor. (Note that 
during an outburst of a transient black-hole system such 
a spherical source may be present, as the central accretor 
may develop a corona.) For a small source at the centre of 



the disc and lying in its plane, the irradiation temperature 
Tirr(i?) is given by 
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(Fukue, 1992). Here r\ is the efficiency of rest-mass energy 
conversion into X-ray heating, j3 is the X-ray albedo, and 
H(R) is the local disc scale height. The minimum accretion 
rate required to keep the disc in the high state is given by 
setting Ti TI (R) = Th, where is Th is the hydrogen ionization 
temperature. Since T always decreases with R, the global 
minimum value Merit i s given by conditions at the outer 
edge i?d of the disc. For the parametrization adopted by 
King, Kolb & Szuszkiewicz (1997), and n — 0.2, this leads 
to 



MZl t (R) = 2.86 x 10 
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where /0.7 is the disc filling fraction / (the ratio of Rd to 
the accretor's Roche lobe) in units of 0.7; mi,mj are the 
accretor and companion star mass in Mq; and 
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Equation (^) is the same as eqn (12) of King, Kolb & 
Szuszkiewicz (1997) apart from the factors /o.7fl, there taken 
as unity. All of the uncertainties over disc thickness, warp- 
ing, albedo etc are lumped into the quantity g. With g ~ 
fo.7 — 1, equation Q appears to be largely successful in 
predicting that systems with reasonably massive (5 — 7Mq) 
black holes and main-sequence companions should be tran- 
sient. By contrast, neutron star systems with main-sequence 
companions should be persistent, as the index of the ratio 
H/R in (Q) is unity, implying more efficient disc irradiation. 
(Equation (^) also implies that lower-mass black hole sys- 
tems might be persistent.) These results suggest that the 
quantity g appearing in (g) cannot be too far from unity. 



3 THE EDDINGTON LIMIT 

Here I concentrate an another aspect of eq. (|^) which does 
not seem to have received much attention. Namely, for large 
enough P, M^.\ t must exceed the Eddington accretion rate 
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M Edd ~ 1 x l(T 8 miM yr" 1 . (5) 

The obvious consequence of eqs. ^) is that for sufficiently 
long orbital periods irradiation will be unable to suppress 
outbursts, as the required central luminosity exceeds the 
Eddington limit, and the system presumably cannot be both 
super-Eddington and persistent. Note that this conclusion 
holds whatever the actual value of the mass transfer rate in 
the particular binary happens to be. Thus we should expect 
to find no persistent LMXBs above a certain critical orbital 
period P cr it ■ For the neutron-star case this was recognised by 
Li & Wang (1998), who found P cr it(NS) ~ 20 d, in agreement 
with observation. For the black-hole case, combining [s]) 
gives 

P crit (BH) * 2.0/ : 7 1 V°- 75 (-^—) ' 75 mY S ml /S d, (6) 

\0.5Ai E dd / 

where we have included a factor (M/0.5M E dd) to allow for 
the fact that the radiation pressure limit for the accretion 
rate M may in practice be below M E< jd- We thus expect 
to find no persistent black-hole LMXBs above this period. 
This is indeed supported by the available data, but hardly 
surprising in view of the difficulty in identifying black holes 
in persistent systems. Note that in high-mass black-hole 
systems such as Cygnus X-l, the powerful UV luminosity of 
the companion star, as well as the small disc size expected 
in a wind-fed system, are both likely to keep the disc hot 
and therefore give a persistent system. 



4 GRO J1655-40 

With g ~ 1 as argued above, the value of P cr it(BH) found 
above is close to the observed period P — 2.62 d of the black- 
hole soft X-ray transient GRO J1655-40 (the nearest peri- 
ods among alternative black-hole systems are P — 6.47 d for 
V404 Cyg and P = 1.23 d for 4U 1543-47). Indeed Kolb et 
al (1997) pointed out the system's proximity to the Edding- 
ton limit during outburst, and Hynes et al. (1998) explicitly 
suggested that no globally steady disc solution might be 
possible for this system with M < M E dd- GRO J1655-40 is 
unusual in at least two respects: 

1. The companion star has spectral type F3 - F6IV and 
mass Mi ~ 2.3M0. On a conventional view, this places it 
in the Hertzsprung gap. The companion star should there- 
fore be expanding on a thermal timescale and thus driving 
a mass transfer rate —M2 ~ 1O _7 M0 yr" 1 (Kolb et al., 
1997). This is well above the appropriate value of M^riti 
making it puzzling that the system is nevertheless tran- 
sient, and far above the mean mass accretion rate of M ^ B = 
1.26 x 10" 10 M Q yr _1 deduced by van Paradijs (1996) from 
observation. Regos, Tout & Wickramasinghe (1998) appeal 
to convective overshooting to increase the main-sequence 
radius of stars of ~ 2Mq . The companion might then be on 
the main sequence rather than in the Hertzsprung gap. This 
implies a slower evolutionary radius expansion, bringing the 
predicted mass transfer rate below Merit- However — M2 is 
still predicted to lie uncomfortably far above M b a - 

2. The system was first detected in an outburst in 1994, 
and had probably been quiescent for at least 30 yr before 
that. Yet two more outbursts followed in the next two years. 

The considerations given here offer explanations for 



both of these unusual features. First, if P > P cr it(BH), the 
system must be transient in some sense, regardless of the ac- 
tual mass transfer rate (cf Hynes et al., 1998). It would there- 
fore not be necessary to appeal to convective overshooting. 
Further, since the system is close to P cr it(BH), it is evidently 
accreting at a value close to the Eddington rate during its 
quasi-steady states (see below), making it natural that M b s 
is much smaller than the predicted mass transfer rate —Mi. 

Second, assuming that the quantity g has a relatively 
constant value close to unity, as argued above, we see from 
([]) that the value of P cr i t (BH) is most sensitive to the filling 
factor / (I consider the effect of dropping the assumption 
g ~ constant below). Thus if / decreases, P cr i t (BH) can in- 
crease above the actual orbital period, allowing irradiation 
to keep the disc in the high state (prolong an outburst) for 
as long as / remains sufficiently small. Hence the unusual 
outburst behaviour of GRO J1655-40 may be explicable in 
terms of the time evolution of the disc size. Encouragingly 
there is some observational evidence (see the discussion in 
Orosz & Bailyn, 1997) that the grazing eclipses seen in the 
optical are time-dependent, just as expected if the disc size 
varies. Moreover Soria, Wu & Hunstead (1999) find evi- 
dence from the rotational velocities of double-peaked emis- 
sion lines that the disc is at some epochs slightly larger than 
its tidal limit. The large resultant torques on the disc suggest 
that this state cannot persist and the disc must eventually 
shrink. 

In fact we do expect / to evolve systematically: in 
the early part of an outburst, the central accretion of low 
angular-momentum material will raise the average disc an- 
gular momentum and thus cause / to increase, hence lower- 
ing Pcrit (BH) and making the system more vulnerable to a 
return to quiescence. However at some stage matter trans- 
ferred from the companion will tend to reduce the angu- 
lar momentum of the outer disc, thus decreasing /, raising 
Pcrit (BH) and allowing irradiation to stabilize the disc in the 
high state. But eventually the disc must grow towards its 
tidal limit, increasing / and thus lowering P cr it(BH) again, 
finally enforcing a return to quiescence. Obviously a full disc 
code is required to follow this sequence in detail and to check 
if it can account qualitatively for the unusual outburst be- 
haviour of GRO J1655-40. 

Clearly, systematic evolution of one or more of the quan- 
tities appearing in g during the outburst could have a similar 
effect in making P cr it(BH) oscillate around the actual orbital 
period P. The most likely alternative candidate is the disc 
aspect ratio H/R, which would appear explicitly with the 
power 1.5 if we substitute for g in ([]). The aspect ratio 
could evolve systematically on a viscous timescale because 
the disc may warp under radiative torques (Pringle, 1996). 
A warp presenting more of the disc surface to the central 
source would tend to stabilize it against a return to qui- 
escence even though the central luminosity was below the 
Eddington limit. Again considerably more work is required 
to check this possibility. 



5 CONCLUSIONS 

I have shown that the Eddington limit implies a critical 
orbital period P cr it(BH) beyond which black-hole LMXBs 
cannot appear as persistent systems. The precise value of 
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Pcrit (BH) is subject to uncertainties expressed by the quan- 
tity g in (^). I have argued that g cannot be very far from 
unity if we are to understand the difference in the stability 
properties of discs in neutron-star and black-hole systems. 
In this case GRO J1655-40 lies much closer to P cr it(BH) 
than any other black-hole system. 

The unusual behaviour of GRO J 1655-40 may result 
from its location very close to P cr it(BH); evolution of the 
disc size or possible radiative warping may move the system 
across the boundary where a sub-Eddington luminosity can 
keep the disc stably in the high state. This system, and 
those at longer orbital periods, probably have central ac- 
cretion rates which are highly super-Eddington during out- 
bursts. Since observed radiative luminosities are mildly sub- 
Eddington, most of this mass must be expelled. Strong sup- 
port for this comes from the observation of P Cygni profiles 
in GRO J1655-40 (Hynes et al., 1998). The superluminal 
jets observed (Hjellming & Rupen, 1995) in an outburst of 
this system may therefore simply represent the most dra- 
matic part of this outflow. 
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